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Abstract 

We outline a procedure for rendering physically-based 
thermal infrared images of simple vegetation scenes . Our 
approach incorporates the biophysical processes that affect 
the temperature distribution of the elements within a scene . 
Computer graphics plays a key role in two respects. First, in 
computing the distribution of scene shaded and sunlit facets 
and, second, in the final image rendering once the tempera- 
tures of all the elements in the scene have been computed. 

We illustrate our approach for a simple com scene where 
the three-dimensional geometry is constructed based on 
measured morphological attributes of the row crop. Statis- 
tical methods are used to construct a representation of the 
scene in agreement with the measured characteristics. 

Our results are quite good. The rendered images ex- 
hibit realistic behavior in directional properties as a func- 
tion of view and sun angle. The root-mean- square error in 
measured versus predicted brightness temperatures for the 
scene was 2. 1 deg C. 


1. Introduction 

Remote sensing requires the interpretation of measured 
sensor response to extract information about earth surface 
features. We can use remote sensing to monitor biogeo- 
physical attributes and assess environmental conditions. 
The interpretation of such signals, however, is a complex 
process dependent on many factors including sensor re- 
sponse function, atmospheric state, and the interplay of the 
probing electromagnetic radiation (passive or active) with 
the underlying scene. 

Models play a strategic role in developing information 
extraction algorithms or in sensor design studies. Tradition- 
ally, these models have been one-dimensional, plane paral- 
lel abstractions. These have proven useful for the current, 
large-scale global monitoring satellite systems such as the 
NOAA Advanced Very High Resolution Radiometer with 


J. R. Ballard, Jr. 
Environmental Laboratory 
US AE Waterways Experiment Station 
Vicksburg, MS 39180 
ballard@orion.wes.army.mil 


its 1 kilometer instantaneous field-of-view. However, with 
the next generation, high spatial resolution satellite systems, 
more attention must be focused on the three-dimensional 
structure of underlying scenes. The new satellite systems 
have spatial resolutions of 1 to 30 meters which is compa- 
rable to many intrinsic scene characteristics, e.g. crop crow 
structure. 

Remote sensing researchers have already begun to draw 
upon traditional computer graphics techniques to aid them 
in their modeling of more complex scenes, particularly in 
the optical reflective regime. Borel, et al. [1] and Goel, 
et al. [3] used the radiosity method developed by Green- 
berg and associates [4] to model the multi-spectral reflec- 
tive response of three-dimensional plant canopies. Smith, 
et al. [9] also employed a statistical version of the radios- 
ity method to compute short-wave absorption of a layered 
canopy for subsequent use in a one-dimensional thermal in- 
frared model. 

While progress has been made in the modeling of 
the infrared response of three-dimensional scenes [6, 7], 
the modeling of the thermal infrared behavior of three- 
dimensional vegetation canopies is still in the beginning 
phases [8]. 

In this paper, we discuss a novel technique for combin- 
ing computer graphics, e.g. ray tracing, with physics-based 
energy budget modeling to predict the thermal infrared ex- 
itance of a plant canopy. We illustrate the method for an 
agricultural scene where meteorological, crop morphologi- 
cal, and thermal infrared measurements were available. 

2. Method 

Figure 1 gives a flow chart of the calculation steps 
required in our approach. First, we generate the three- 
dimensional geometry of the scene using simple measure- 
ments of the canopy including crop row spacing, height, leaf 
slope distribution, and a measure of the density of leaf area 
per unit ground area, the Leaf Area Index (LAI). The details 
are described below. 



We then use ray tracing to determine all of the elements 
in the scene which are either directly illuminated, i.e., sun- 
lit, or in shade. For each of these sunlit or shaded elements, 
we separately run a thermal model to calculate their energy 
budget and to solve for the corresponding leaf or soil tem- 
peratures. We subsequently, again, use ray tracing to render 
the thermal infrared scene by tracing the thermal infrared 
flux, emitted by each element at its estimated temperature, 
which is actually received by a sensor at a specified view 
direction. The thermal model is driven by local meteoro- 
logical conditions including local wind speed, relative hu- 
midity, air temperature, and short-wave illumination and re- 
quires only simple material properties. 

We used a nested approach to simulate the geometry 
of the scene based on measured canopy parameters. We 
modeled a sample plot of 19 m by 19 m for which canopy 
measurements, meteorology, and thermal infrared measure- 
ments were available. We first generated twenty-five unique 
19 m rows with a canopy width of 0.46 m, height of 0.8 m, 
total row spacing of 0.76 m including foliage and gap be- 
tween rows, and a single sided leaf area index of 3.0 cor- 
responding to the measured row spacing, height, and LAI. 
A spherical leaf slope distribution and uniform leaf azimuth 
distribution were used to assign orientations to the individ- 
ual leaf facets. Our primitive leaf facet size was 5 cm by 5 
cm providing a leaf area of 0.0025 m 2 . Thus, 10,488 leaf 
facets were required to match the measured LAI of 3.0 dis- 
tributed over the 19 m by 0.46 m row area. The twenty-five 
rows were instanced together at a row spacing of 0.3 m to 
simulate a corn foliage plot measuring 19 m by 19 m con- 
taining approximately 262,200 facets. 

Figure 2 is a visual simulation of a hemispherical picture 
of the scene. The sun was 73.2 deg above the horizon and 
solar azimuth was - 8.8 deg measured from the south. The 
sunlit and shaded components of both canopy and soil el- 
ements are apparent. From such images, we can compute 
a shade mask for all elements within the scene in order to 
apply the appropriate energy balance calculations discussed 
below. 

Individual leaf energy budgets are calculated as follows 
in terms of net radiation, R n , sensible heat, H, and latent 
heat, LE, i.e. 

R n = LE + H (1) 

A standard resistance formulation [2] is used for the calcu- 
lation of latent and sensible heat in terms of the unknown 
leaf temperature, T vy i.e. 
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where T a is the air temperature, h veg is the relative humid- 
ity of the vegetation air space, e* is the saturation water 


vapor pressure at canopy temperature T Vi e a the vapor pres- 
sure of air at air temperature T a , p is the air density, c p 
the air specific heat, and 7 the psychromatic constant. Fi- 
nally, r& is the canopy boundary layer resistance and r c is 
the bulk stomatal resistance. Canopy stomatal resistance 
and the boundary layer resistance depend upon meteorolog- 
ical driving variables. 

Net radiation to a surface, R n , is computed as the sum of 
the absorbed long-wave and short-wave radiation minus the 
outgoing long-wave emitted thermal infrared flux. Separate 
calculations are required for sunlit versus shaded surfaces. 
For sunlit leaves, the absorbed shortwave flux is highly de- 
pendent on leaf inclination angle or the angle between the 
leaf and the incident solar direction. Long-wave emission is 
a function of the unknown leaf temperature and is given by 
the Stefan-Boltzmann relation. 

Thus, a system of non-linear energy budget equations is 
solved for individual leaf temperatures for all possible leaf 
orientations corresponding to varying exposure of the indi- 
vidual leaves to direct solar flux. For sunlit leaves, one de- 
gree increments are used for leaf inclination classes. For the 
case simulated, sunlit leaf temperatures temperatures range 
from 24 to 32 deg C with warmer temperatures occurring 
more frequently. For shaded leaves, the average leaf tem- 
perature was 26.5 deg C. Average soil sunlit and shaded 
temperatures corresponding to the measured meteorologi- 
cal conditions were 55 and 46 deg C respectively. 

For the ray-tracing and rendering calculations we em- 
ployed the RADIANCE program developed by Ward [10]. 
Small modifications were required to handle the emission 
from both sides of the leaf facets. 

3. Results 

Figure 3 displays a simulated thermal infrared corn im- 
age corresponding to Figure 2 rendered at one-quarter me- 
ter pixel size. High LAI or dense vegetation yields cooler 
temperatures in the thermal image compared to the warmer 
soil facets which are primarily visible looking straight down 
into the image. 

The brightness temperature of a surface is defined as the 
equivalent black-body temperature of the surface, i.e. a sur- 
face with emissivity 1 . 0 , that would yield the same total 
outgoing long-wave flux as the measured surface. Ground 
measurements of brightness temperature and meteorologi- 
cal driving variables were available for the 19 m by 19 m 
sample area. As a check on our energy budget formula- 
tions, we calculated net radiation for this area using Equa- 
tions ( 1 ) and ( 2 ) and simulated the three-dimensional ther- 
mal infrared response for the test plot corresponding to 600, 
800, 1000, 1200, 1400, 1600, 1800, and 2000 hours. Fig- 
ure 4 shows a comparison of the measured versus predicted 
brightness temperatures for these 8 time periods as calcu- 



lated from the hemispherically averaged outgoing thermal 
infrared flux, assuming an emissivity of 1.0. The root- 
mean-square error was 2.1 degrees C with an r 2 of 0.93. 
The error bars correspond to the estimated confidence in- 
tervals from the regression line fit to the predicted versus 
measured values. 

4. Summary and Conclusions 

In this paper we have given a general outline for ren- 
dering thermal infrared images of simple vegetation scenes 
based on physical principles. Our approach incorporates 
the biophysical processes that affect the temperature distri- 
bution of the vegetation elements within a scene. The ther- 
mal model is driven by readily available, standard meteo- 
rological measurements including local wind speed, rela- 
tive humidity, air temperature, and short-wave illumination. 
Computer graphics plays a key role in two respects. First, 
in computing the distribution of shaded and sunlit facets 
within the scene and, second, in the final image rendering 
once the temperatures of all the elements in the scene have 
been computed. 

We illustrated our approach for a simple corn scene 
where the three-dimensional geometry was constructed 
based on measured macroscopic morphological attributes of 
the row crop— viz crop height, row spacing, leaf-area-index, 
etc. Statistical methods were then used to randomly assign 
individual leaf facets with the proper orientations within 
rows and to construct a representation of the scene in agree- 
ment with the measured characteristics. 

Our results are quite good. The rendered images ex- 
hibited realistic behavior in the directional properties as a 
function of view and sun angle [8]. The root-mean-square 
error in measured versus predicted brightness temperatures 
for the scene was 2. 1 deg C. 

Nevertheless, there are several limitations in our method. 
First, our vegetation scene only consisted of leaves with no 
woody material, e.g. trunks. The time constant for indi- 
vidual leaves to reach thermal equilibrium is less than 30 
seconds and our equilibrium energy-budget approach was 
thus warranted. In soil, we used a simple trick based on re- 
sults reported in the literature, i.e. modeling conduction as a 
fraction of net radiation (possibly with time lag). In forests, 
however, the woody material, e.g. trunks and branches, 
must be accounted for. It is more challenging to incorpo- 
rate the geometry of these elements. We are experimenting 
with Lindenmayer systems [5]. 

Finally, we imposed a static wind field on the scene, con- 
sistent with average wind speed measurements and a loga- 
rithmic height profile. More realistic modeling of high fre- 
quency phenomena, e.g. leaf flutter, opens up the pandora’s 
box of turbulence, a subject in its own right. We are in- 
vestigating all of these effects and further work is clearly 


warranted. 
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Figure 1. Flow chart of model calculations 



Figure 2. Simulated hemispherical visual image 
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Figure 3. Simulated hemispherical thermal infared image 
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Figure 4. Comparison of predicted versus measured canopy brightness temperature 




